Where the granular flows bend 
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ABSTRACT 

Based on IMaX/SuNRlSE data, we report on a previously undetected phe- 
nomenon in solar granulation. We show that in a very narrow region separating 
granules and intergranular lanes the spectral line width of the Fel 5250.2 A line 
becomes extremely small. We offer an explanation of this observation with the 
help of magneto-convection simulations. These regions with extremely small line 
widths correspond to the places where the granular flows bend from mainly up- 
flow in granules to downflow in intergranular lanes. We show that the resolution 
and image stability achieved by IMaX/SuNRlSE are important requisites to de- 
tect this interesting phenomenon. 

Subject headings: Sun: granulation - Sun: observations - Sun: spectral line 
formation 



Introduction 



Solar granulation is one of the most studied phe nomena in the lo wer solar atmosphere 
and its properties are thought to be well understood (IBray et al.lll984j ). Spectral and imag- 
ing observational data have allowed solar granulation to be studied in detail. Commonly 
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considered parameters include sizes and lifetimes of granules, velocity and intensity flue- 
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). At the same time, state of the art simula- 



tions of solar convection and magneto-convection reproduce very satisfactorily many obser- 
vational granulation properties, which has made it possible to explain self-consistently the 
mean shapes of spectral line profiles (INordlund et al.l 120091 ). Simulations have shown that 
parameters as micro- and macroturbulent velocities (used widely in the past to match the 
spectral lines shapes to observations) are not necessary any more as the spectral lines are 
broadened by the velocity field of the convection motions. 

In this Letter, we report on spatially resolved observations of the line broadening by 
the granular velocity field (or, equivalently, the spatially resolved action of solar micro- 
turbulence). The data used in the Letter were ob tained with the Imaging Magnetograph 



experiment (IMaX; see lMartmez Pille 



telescope in June, 2009 (IBarthol et al. 
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2010 
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Solanki et al.l 1201 01 ). We study the spatial vari- 



ations of the full width at half maximum (FWHM) of the Fel 5250.2 A line. We found 
that the FWHM is larger in intergranular lanes and smaller in granules. But the most 
interesting finding is that the FWHM is particularly small in narrow lanes located at gran- 
ular/intergranular borders. We fou nd only few work s that discuss the relatio n between the 



line FWHM and granulation (e.g. iNesis et al.l Il992t lHanslmeier et al.l 120081) . pointing out 



that the FWHM is anti-correlated with brightness. The observations of INesis et al . (1992) 
were interpreted in terms of supersonic flows, which was questioned by lSolanki et al.l ( 119961 ) . 
We did not find any previous mention of the detection of the narrow regions with extremely 
small FWHM adjacent to intergranular borders. The aim of this Letter is to describe their 
properties and to offer an explanation for their origin. 



2. Observed line width 



The IMaX instrument and the d etails of the data acquisition and reduction are de- 
scribed in lMartmez Pillet et al.l (120101 ). Here we use a time series of images of about 23 min 
duration taken near the solar disc center (UT 00:37). The observational sequence consists 
of five polarized filtergrams taken at ±40, ±80 and ±227 mA from the Fel 5250.2 A line 
center. We use Level 2 data (reconstructed by deconvolution) . The spatial resolution of 
the reconstructed data is estimated to be 0.15. The field of view is about 33" x 33" (after 
removing the apodized areas), with a pixel size of 0.055. We use the following parameters 
extracted from the Stokes / data: continuum intensity, line of sight velocity, and FWHM. 
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Fig. 1. — (a) Image of the full width at half maximum (FWHM) of the Fel 5250.2 A line 
in the IMaX field of view; the grey scaling ranges from 90 to 150 mA; the mean FWHM is 
equal to 115 mA; (b) same image after application of the unsharp masking filter. Panels (c) 
and (d) give an enlarged view of the area of 6 Mm x 6 Mm marked by a white square in 
panels (a) and (b). Blue contours in panel (d) correspond to normalized continuum intensity 
Ic/Ic — 1- Red dashed contours mark locations with zero line of sight velocity. 
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Fig. 2. — Left panel: temporal variations of the FWHM in the IMaX data along a virtual 
slit placed at one selected horizontal location; the grey scaling ranges from 90 to 150 mA; 
the mean FWHM is equal to 115 mA. Middle panel gives the same image after application of 
the unsharp masking filter. Right panel: histograms of the LOS velocity in the IMaX data at 
the locations where the continuum intensity is above unity, i.e. in granules (red solid line); 
at the locations where the continuum intensity is below unity, i.e. in intergranular lanes (red 
dashed line); at the locations with enhanced FWHM marked white in Fig. [lb (black dashed 
line); and at locations with extremely small FWHM marked black in Fig. [Tb (black solid 
line). 
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The latter two parameters were calculated from Gaussian fits to the five intensity wavelength 
points. The statistical errors of the fitting procedure were within the range of 20-40 m s _1 
for the velocities and a few mA to 10 mA (at most) for the FWHM. 

Figured] shows an image of the FWHM at one moment of time. This image demonstrates 
clearly that the width of the Fel 5250.2 A intensity profile follows the granulation structure. 
The FWHM increases in intergranular lanes (appearing as bright regions) and becomes 
smaller in granules. A more careful inspection of the image reveals that adjacent to almost 
each intergranular lane there are narrow black regions (that also have elongated lane-like 
shapes) of smaller widths. The average FWHM of the Fel 5250.2 line in the granular regions 
is of 120 mA, in intergranular regions is of 140 mA, but only 100 mA at the borders between 
granules and intergranular lanes. The unsharp mask filteio applied to this image enhances 
the gradients and makes the locations with extremely small FWHM more evident (Fig. [Ib). 
Panels c and d with enlarged view of the region of 6 Mm x 6 Mm size (marked by the white 
squares in Figs. [T^i and b) give some very clear examples of the dark lanes surrounding bright 
intergranular spaces. These dark features sometimes appear only on one side of the bright 
intergranular lane, and sometimes on both sides. In many cases they appear in regions where 
the normalized continuum intensity J/J c ~ 1, i.e. just at the borders between granules and 
intergranular lanes. The width of these dark features is about 3 — 4 IMaX pixels, i.e. ~o''2. 
Their lifetime is the same as the lifetime of the intergranular lane they are attached to, see 
Fig. |2] (left and middle panels) . These dark lanes in the FWHM image appear also in the 
non-reconstructed data (with less contrast), proving that they are not artifacts of the phase 
diversity reconstruction procedure. 



3. Origin of the dark lanes 



To find the origin of these dark lanes in FWHM images, we have performed spectral 
l ine synthesis of th e Fel 5250.2 A line in a snapshot of magneto-convection simulations 
( IVogler et al.ll2005f ). We considered a simulation run with a bi-polar magnetic field structure 
and an average unsigned magnetic field in the box of 30 G, corresponding to a quiet solar 
region (for more details on this simulation snapshot see Khomenko et al.l 120051 ). The grid 
size of the simulations is about 20 km The synthesis was done in LTE using the SIR code 
( jRuiz Cobo fc Del Toro Iniestal 119921 ). No additional line broadening mechanisms, such as 
macro- or microturbulence were introduced. 



We found that the areas with extremely narrow intensity profiles adjacent to intergran- 



The algorithm that works by enhancing the contrast between neighboring pixels in an image 
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Fig. 3. — (a) Image of the FWHM of the Fel 5250.2 A line obtained from simulations 
degraded to IMaX resolution, the grey scaling ranges from 90 to 150 mA; the mean FWHM 
is equal to 112 mA; (b) same image after application of the unsharp masking filter; blue 
contours correspond to l c jl c = 1; red dashed contours correspond to zero line of sight 
velocity; green lines indicate locations considered in Fig. 5; (c) image of the FWHM of the 
Fel 6301.5 A line obtained from simulations degraded to Hinode resolution; the grey scaling 
ranges from 123 to 160 mA; the mean FWHM is equal to 132 mA. 
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Fig. 4. — Left panel: image of the FWHM obtained from simulations when the spectral 
synthesis is done after setting horizontal variations of thermodynamics variables to zero; the 
grey scaling ranges from 90 to 150 mA; the average value is 112 mA. Right panel: the same 
for the case that the synthesis is performed after setting velocities to zero in the simulation 
snapshot; the grey scaling ranges from 98 to 105 mA, the average value is 100 mA. 
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ular lanes also appear in the simulations. For an adequate comparison with the observations 
we have degraded the simulation spectra in the following way. Firstly, we performed a convo- 
lution of each monochromatic image with an Airy function imitating the action of a 100 cm 
telescope. Secondly, we convolved the spectra in wavelength with the transmission profiles 
of the IMaX etalon. Finally, we fit a Gaussian function to the five wavelength points and 
calculated the width of the Gaussian profiles as a proxy to the intensity profile width, as 
well as velocity, similar to the observations. The resulting image of the simulated FWHM is 
shown in Figure^. 

The FWHM image shows structures similar to the observed ones, i.e. larger FWHM 
in intergranular lanes and very low FWHM in the narrow regions adjacent to lanes. Note 
that the average value and the range of variations agree very well in the simulations and 
observations. Similar to the observations, the dark lanes in the FWHM image correspond 
to the normalized continuum intensity around unity, i.e. they mark the transition between 
granules and intergranular lanes (Fig. [3b). Their width is a few pixels (~o'.2). 

Due to the small width of these dark lanes, about o'.'2, high resolution and image stability 
are required for their detection. Figure [3b shows the FWHM after the simulated spectra were 
degraded to Hinode resolution and pixel size (note that the Fel 6301.5 A is used in this case). 
This image preserves brightenings in intergranular lanes corresponding to large FWHM, but 
the adjacent dark lanes do not show up at this resolution. On the one hand, this comparison 
demonstrates that IMaX/SuNRlSE data quality is necessary to detect these features. On 
the other hand, it pr oves that the instrument reached indeed very high spatial resolution 



flBerkefeld et al.ll2010j ). 



The IMaX spectral line of Fel 5250.2 at A is sensitive to the velocity broadening, but 
also to magnetic fields via the Zeeman effect, and to temperature. To find out which effect 
is primarily responsible for the spatial variations of the FWHM, we performed additional 
calculations. We repeated the spectral synthesis for three cases: (1) magnetic field was set 
to zero, preserving the variations in the rest of the variables; (2) horizontal variations of the 
temperature and gas pressure were set to zero using for all the points the mean stratification, 
but preserving variations of the velocity and magnetic field; (3) velocity was set to zero at 
all points, while preserving variations in the rest of the variables. 

We found that the magnetic field's influence on the horizontal variations of the FWHM 
is negligible except for a few isolated locations in intergranular lanes where the field reaches 
kG values. There the FWHM is additionally enhanced by the Zeeman broadening. Figure 
H] gives the variations of the FWHM in the cases (2) and (3). Case (2), where only velocity 
variations are preserved, looks almost identical to the original case in Fig. [3K While case 
(3), where only variations of the thermodynamic variables are preserved, looks very different, 
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showing almost an opposite behaviour to the case (2). Note that, curiously, in the case (3) 
the broadening due to magnetic field becomes evident in some intergranular lanes where the 
field is strongest (locations at X = 1 - 2 Mm; Y = 1 - 2 Mm and X = 4 Mm; Y = 3 Mm 
in the right panel of FigHJ). The broadening due to magnetic field is much weaker compared 
to the velocity broadening, but is comparable in magnitude to the thermal broadening. The 
average value of the FWHM and its standard deviation is much lower for the case (3) as 
compared to the case (2) and to the observations. Thus, we conclude that velocity broadening 
is the prime responsible for the horizontal variations of the FWHM. 

Figure [5] shows several examples of the velocity field at the locations marked by letters in 
Fig. [3b. At these locations the red line crosses an intergranular lane with enhanced FWHM 
as well as the transition between granule and intergranule with very low FWHM. As follows 
from Fig. [51 all these locations have features in common. At the center of an intergranular 
lane, in the deep photospheric layers, the vertical velocity is mainly a downflow with a 
magnitude reaching about 3 — 4 km s" 1 . The intergranular lanes are also characterized by 
a strong vertical velocity gradient. In a few hundred kilometers, the velocity decreases from 
3 — 4 km s _1 to zero, following the change of its sign at the top of the atmosphere. These 
strong velocity gradients have immediate i mplications on the width of the Fel 5250.2 A line 



profiles, whose core forms around 400 km (IGurtovenko &: Kostikj|l989f ). The action of these 



vertical velocity gradients over the height range of formation of the line is similar to the action 
of microturbulence velocity, producing additional broadening of the profiles. We conclude 
that the large FWHM in intergranular lanes is produced by a microturbulence-like effect 
due to stro ng vertical v e locity gradients, unlike the mechanism based on super-sonic flows 



proposed by lNesis et al.l (119921 ) . Similar conclusion, as for the enhancement of the FWHM in 



intergranular lanes, was previously reached in the work of iGadun et al.l (119971 ) from studying 



the line-parameter variations in the two-dimensional hydrodynamic simulations. 

Figure [5] shows that the regions with extremely low FWHM correspond to the locations 
where the granular flow bends. These locations are the transitions from granules with up- 
flowing velocities to intergranular lanes with downflows. As the velocity vector bends there, 
the vertical velocity component is close to zero and the vertical velocity gradient is very 
small as well. It makes the velocity broadening at these transition locations extremely low, 
thus producing narrow profiles, as detected in the IMaX observations. In Fig. [3b the red 
dashed contours mark the location with zero vertical velocity in simulations. It shows that 
almost always these locations coincide with dark lanes in FWHM image. 

The IMaX data only have five points in wavelength, making it impossible to check the 
conclusions about the velocity gradients, obtained from simulations. However, the absolute 
values of the LOS velocities are reliably obtained from these data. According to simulations, 
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the regions with extremely low FWHM correspond to locations with almost zero vertical ve- 
locities. In Fig. [2] (right panel), we plot the LOS velocity histograms obtained from the IMaX 
data over the whole time series. This figure shows the existence of the three distinct families 
of features. Granules, selected as locations with I c /I c > 1, have the velocity histogram peak 
at —0.6 km s _1 . Intergranular lanes, selected as locations with I c /I c < 1, show downflowing 
velocities with the histogram peak at 0.4 km s -1 . This histogram almost coincides with 
the one obtained by selecting locations with enhanced FWHM given by white regions in 
Fig. [Tb. Finally, there is a third distribution corresponding to the borders between granules 
and intergranular lanes. It is obtained by selecting locations with low FWHM given by the 
black regions in Fig. [lb. Its maximum appears at zero LOS velocity. This observational 
confirmation supports the conclusion based on simulations that the regions with extremely 
small FWHM correspond to locations where the velocity vector bends, producing a zero LOS 
velocity component. 

4. Conclusions 

In this Letter we have investigated the spatial distribution of the full width at half max- 
imum of the Fe I 5250.2 A intensity profiles in the IMaX/SuNRlSE data. We have shown that 
the FWHM becomes large in intergranular lanes and is very small in narrow o'.'2 wide strips 
adjacent to the lanes. The lifetime of these narrow regions is the same as the intergranular 
lanes they are attached to and they evolve in parallel. Based on the similarity between 
IMaX observations and simulations of magneto-convection, we have offered an explanation 
for the formation of these profiles with small FWHM. We suggest that these profiles appear 
in regions where the granular flow bends. There, the velocity broadening is extremely small, 
as the LOS velocity component is close to zero and has very small gradients. This conclusion 
is supported by the analysis of the LOS velocity histograms in the IMaX data. We claim 
that the SuNRlSE/IMaX resolution and image stability is the necessary condition to detect 
these narrow transition zones between granules and intergranular lanes. 
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Fig. 5. — Examples of the velocity field in the simulations at locations marked (A), (B) 
and (C) in Fig. [3b. The color coding indicates vertical velocity, with red being a downflow 
and blue being an upflow, both saturated at 3 km s -1 . The white and black dotted lines 
indicate the effective formation height of the continuum and line core of the Fel 5250.2 A, 
respectively, estimated from the calculation of the Response Functions. The white solid line 
shows horizontal variations of FWHM in arbitrary units. 



